The Anderson transition: time reversal symmetry and universality 
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We report a finite size scaling study of the Anderson tran- 
sition. Different scaling functions and different values for the 
critical exponent have been found, consistent with the exis- 
tence of the orthogonal and unitary universality classes which 
occur in the field theory description of the transition. The 
critical conductance distribution at the Anderson transition 
has also been investigated and different distributions for the 
orthogonal and unitary classes obtained. 



It is now widely accepted that the metal- insulator transi- 
tion for non- interacting electrons, the so called Anderson 
transition Q, is a continuous phase transition in which 
static disorder plays a role analogous to temperature in 
thermal phase transitions. The field theoretical formu- 
lation of the problem 0, though not making reliable 
predictions for the critical exponent does indicate 

that it should be possible to describe the critical behavior 
within a framework of three universality classes: orthog- 
onal, unitary and symplectic. However, recent work put 
this idea in question. It was found that the scaling func- 
tions for systems with orthogonal and unitary symmetry 
could not be distinguished at an accuracy of a few percent 
H . Nor could the values of the critical exponent for the 
three universality classes be reliably distinguished 
This, together with recent work on the statistics of en- 
ergy levels in the vicinity of transition, has prompted the 
suggestion Jl(| that the universality classes predicted by 
the field theory do not correctly describe the Anderson 
transition. 

The two important symmetries in the field theory are 
time reversal symmetry (TRS) and spin rotation sym- 
metry (SRS). The system is said to be in the orthogonal 
universality class if it has both SRS and TRS, in the uni- 
tary class if TRS is broken and in the symplectic class 
if the system has TRS but SRS is broken. The relevant 
terms in the Hamiltonian are a coupling to an applied 
magnetic field, which breaks TRS, and the spin orbit in- 
teraction, which breaks SRS. 

Here we focus on the breaking of TRS by a constant 
applied magnetic field. We report the results of Monte 
Carlo studies which indicate that the critical behavior, 
at least as far as orthogonal and unitary symmetries are 
concerned, is in accord with the conventional universal- 
ity classes. By carrying out a precise study of the fi- 



nite size scaling of the electron localization length, which 
is analogous to the correlation length in thermal transi- 
tions, we have clearly differentiated the scaling functions 
for the orthogonal and unitary universality classes. We 
have also estimated the critical exponents and calculated 
confidence intervals for these estimates. We find a sta- 
tistically significant difference of about 12% between the 
values of the critical exponent in the orthogonal and uni- 
tary classes. 

To reinforce the above conclusion we have also sim- 
ulated the critical conductance distribution of a disor- 
dered mesoscopic conductor. Following the discovery of 
universal conductance fluctuations it was realized that 
the conductance of a phase coherent system is not self- 
averaging. Extrapolating from the metallic regime, it 
seems that the conductance fluctuations at the critical 
point should be of the same order of magnitude as the 
mean conductance and that the full conductance distri- 
bution should be a more useful characteristic of the criti- 
cal point. As we approach the critical point the localiza- 
tion length diverges and the system becomes effectively 
self similar, we then expect the conductance distribution 
to become independent of system size and depend only on 
the universality class. This expectation was borne out in 
our study where we found different critical conductance 
distributions depending on whether or not TRS is bro- 
ken. 

The model Hamiltonian used in this study describes 
non- interacting electrons on a simple cubic lattice. With 
nearest neighbor interactions only we have 

<f]H\r> = V(f), 

< r\H\f-x> = 1, m 
<r\H\r-y> = 1, 1 > 

< r\H\f — z > = exp(— i1~K<f>x), 

where x, y and z are the basis vectors of the lattice. The 
electrons are subject to an external magnetic field ap- 
plied in the y direction whose strength is parameterized 
by the flux (f>, measured in units of the flux quantum 
h/e, threading a lattice cell. The on site energies of the 
electrons {^(r)} are assumed to be independently and 
identically distributed with probability p(V)dV where 

p(V) = 1/W \V\ < W/2, 
= otherwise. 

The critical point, scaling function and the value of the 
critical exponent for several values of <p and Fermi energy 



1 



E f (see Table J) are determined by examining the finite 
size scaling [11] of the localization length A for electrons 
on a quasi- Id dimensional bar of cross section Lx L. The 
localization length A = X(Ef, <fi, W, L), defined by 



A" 



lim 



< - In g(L z ) > 
2L, 



(2) 



where L z is the length of the bar and g(L z ) is the con- 
ductance of the bar measured in units of (e 2 /h), can be 
evaluated by rewriting the Schroedinger equation as a 
product of transfer matrices; A -1 can then be determined 
to within a specified accuracy using a standard technique 
p2[ . The accuracy used here ranges between 0.1% and 
0.2%. 

On intuitive grounds it has been argued [|l3| that we 
should observe orthogonal scaling when L <C Lh and 
unitary scaling when L ^> Lh where Lh is the magnetic 
length. For the lattice model (|l|) Lh = ^l/2w<p. For 
the smallest system size L = 6 used here this criterion 
yields a crossover flux C ~ 1/226. We thus expect to see 
clear unitary scaling behavior for cases Ua and Ub listed 
in Table § 
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FIG. 1. The scaling of the data for the parameter sets listed 
in Table Ijj. The lines are the scaling functions given by (^) 
and (|^). Different scaling functions for the orthogonal and 
unitary universality classes can be clearly distinguished. 

When the dimensionless quantity A = X/L is plot- 
ted against disorder W for different cross sections L the 
curves are found to have a common point of intersec- 
tion; this indicates the occurrence of the metal- insulator 
transition at a critical disorder W c in the 3c? system which 
would be obtained by letting L — > oo. Detailed analysis 
of the data is based on the following assumptions: first 
that for L finite A is a smooth function of W and L, 
second that the data obey a one parameter scaling law 



where w = (W c — W)/W c and the subscript ± refers to 
w > and w < 0, and third that the length £ appearing 
in the scaling law has a power law divergence close to 



W c of the form £(w) 



w 



This relation defines 



the critical exponent v and introduces two arbitrary con- 
stants £ + (w > 0) and £~ (w < 0). According to the 
Wegner scaling law |Q v is related to the critical expo- 
nent s associated with the conductivity by s = (d — 2)i> 
so that v = s in d = 3. The above assumptions imply 
that it should be possible to fit the data to 



Afi t (W, L) = A c + J2 A n L n '»w r ' 



(4) 



In practice we have truncated this series at n = 3. The 
relation between (||) and ([|) can be made apparent by 
writing 



/+=A c + E n=1 a«(i/0 



An (? 



(5) 



/- = a c + e„=i b n {L/er /v , b n = (-i)m„ (n n/ " 



(6) 



In principle £ + and £~ should depend on energy and flux 
though the "amplitude ratio" £ + may be universal 
|p^| so that £ + and £~ may not be independent. Their 
absolute values cannot be determined using the present 
method. No relative variation as a function of energy 
and flux, which would be apparent in the simulation, was 
detected. Therefore for convenience we set £ + = £~ = 1. 
The most likely fit is determined by minimizing the \ 2 ~ 
statistic 



E 



A, 



A Rt (W m ,L m ) 



(7) 



where the summation m is over all data points and a m 
is the error (standard deviation) in the determination of 
the mth data point. After being fitted the data are re- 
plotted against L/£ to check that they obey the scaling 
law @. 

We also need to determine the goodness of fit Q and 
confidence intervals for the fitting parameters. The good- 
ness of fit measures the credibility of the fit; Q > 0.001 
is often regarded as acceptable in other applications [|l6| . 
We have checked that the numerical procedure used to 
estimate the localization lengths does so with an error 
which is approximately normally distributed. If we ig- 
nore the presence of any systematic corrections to scal- 
ing in the data, this permits the use of the x 2 likely-hood 
function to determine the "best fit" and the estimation 
of Q from the \ 2 distribution with M — N degrees of 
freedom. 



A(W,L) = /±(L/£H), 



(3) 
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FIG. 2. The critical conductance distributions for parame- 
ter set Oa listed in Table 

The confidence intervals for the fitted parameters were 
estimated in two independent ways: first from the Hes- 
sian matrix obtained in the least squares fitting proce- 
dure and second using the Bootstrap procedure described 
in Jl6| . In the latter method the original data are repeat- 
edly randomly sampled (with replacement) and fitted. 
This provides an independent check of the distribution 
of the fitted parameters. Both methods gave approxi- 
mately the same results. We chose to present the errors 
as 95.4% marginal confidence intervals as given by the 
Boostrap method. 

The results are summarized in Tables | and 0. Given 
the confidence intervals the probability that the values 
of the critical exponent for the orthogonal and unitary 
universality classes are the same seems to be negligible. 
Different values of A c can also be distinguished confirm- 
ing what is clearly evident in Fig. [l] that the orthogonal 
and unitary data scale differently. We conclude that the 
scaling function is sensitive to the breaking of time re- 
versal symmetry. 

We now turn to the conductance distribution. We con- 
sider a cubic "sample" of side L attached to semi- infi- 
nite leads on two opposite faces. The disorder W is set 
to zero in the leads, the Fermi energy and magnetic field 
are constant throughout. The zero temperature linear 
conductance G — (e 2 /h)g can be obtained from g = trtt^ 
where t is the transmission matrix of the sample. The t 
matrix can be related to a Green's function [17 18 which 
can be determined iteratively. 

For each set of parameters we have calculated the con- 
ductances for an ensemble of 100, 000 realizations of the 
random potential. Some typical results are presented in 
Fig. |[ Our analysis of these results is based on the 
generalization of (||) p(g) — p{g,L/£). At the critical 
point £ diverges and we should obtain a universal critical 



conductance distribution p c (g) which is independent of 
L. The scale invariance of p c (g) is clearly demonstrated, 
at least for the range of system sizes studied, in Fig. 
^. We found a similar scale invariance for all the cases 
listed in Table |[ In Fig. || we have plotted the critical 
conductance distributions obtained and in Table III tab- 



ulated some averages of these distributions. The results 
are consistent with the existence of distinct orthogonal 
and unitary critical conductance distributions. 

We now discuss the general features of p c (g) focusing 
on the orthogonal universality class and making a com- 
parison with the critical distribution obtained in the e 
expansion in the field theory. The nth cumulant c n of 
p c (g) for a d dimensional cube of linear dimension L is 
M 



(L/iy 



n 
n 



(8) 



where e = d — 2, no is an integer of order 1/e and I 
is the elastic mean free path. As described in j2(J] it is 
possible, under certain assumptions, to derive p c {g) from 
this. Extrapolating to d — 3 we find 



1 



Pe(g) = S{g) + 



7T 

2e 



U( 2 )+ 2! [M * U]( 2 ) + 



(9) 



where u(x) = 4a; -3 exp(— 2/x), * denotes the convolu- 
tion and e=2.71828... The series (||) is easily handled 
numerically. The result is shown in Fig. |^. 
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FIG. 3. The distribution of g at the critical point. Orthog- 
onal and unitary distributions can be clearly distinguished. 
The critical conductance distribution obtained in the e ex- 
pansion is also shown. 

The most obvious feature is that p c (g) is not peaked 
about its mean value < g >. The conductance fluctua- 
tions, as measured by the standard deviation a g , are of 



3 



the same order of magnitude as < g >. If we compare 
with p e (g) we find that we have a good approximation 
to the central region of the distribution function but a 
rather bad approximation to its tails. The large g tail 
of Pe(g) decays as 1/g 3 which means that all cumulants 
higher than c\ diverge. This is reflected in (Q) where 
these cumulants are not universal but depend on I and 
L. We could find no evidence of this behavior, however, 
in the simulation; as seen in Fig. || there is sharp decay 
of p c (g) above g — 2 and the higher cumulants, at least 
as far as n = 4, seem to have universal values. 

Another way to look at the critical distribution is to 
change variables to In g (see Fig. [|). While the distri- 
bution is certainly not of Gaussian form, it does show a 
central tendency. 
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FIG. 4. The critical distribution after changing variables. 

In conclusion we have presented numerical evidence 
which, we think, confirms that the critical behavior at the 
Anderson transition is sensitive to perturbations which 
break time reversal symmetry. 
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Ef 


<t> 


L 


[Wmin, W ma x] 


M 


N 


Xmin 


Q 


Oa 








{6,8, 10,12} 


[15, 18] 


86 


6 


79 


0.50 


Ob 


0.5 





{6,8, 10} 


[15.6, 17.6] 


73 


6 


49 


0.96 


Ua 





1/3 


{6,9, 12} 


[17,20] 


55 


6 


35 


0.94 


Ub 





1/4 


{8, 10, 12, 14} 


[17.5, 19.5] 


66 


6 


65 


0.31 



TABLE I. The parameters used in the numerical study: 
Fermi energy Ef, flux </>, system sizes L, the range of dis- 
order, the number of data points M, the number of fitting 
parameters N , the value of x 2 for the best fit Xmm an d the 
goodness of fit Q as estimated from the \ 2 distribution. 



Label 


W c 


A c 


V 


Oa 
Ob 
Ua 
Ub 


16.448 ± .014 
16.442 ± .018 
18.316 ±.015 
18.375 ±.017 


0.5857 ± .0012 
0.5862 ± .0015 
0.5683 ± .0013 
0.5662 ± .0016 


1.59 ± .03 

1.60 ±.06 
1.43 ± .04 
1.43 ± .06 



TABLE II. The best fit values of the critical parameters 
together with their 95.4% confidence intervals as estimated 
using the bootstrap method. 
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Label 


<5 > 


°" 9 


< ln S > 


fin g 


OaL = 


10 


0.579 


0.598 


-1.206 


1.319 


Ob L = 


10 


0.578 


0.598 


-1.207 


1.318 


UaL = 


12 


0.395 


0.499 


-1.869 


1.590 


Ub L = 


10 


0.395 


0.499 


-1.864 


1.583 



TABLE III. The means and standard deviations of the crit- 
ical distribution of g and In g. 



